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ABSTRACT 

■ We study the role of rotating sunspots in relation to the evolution of vari- 

ous physical parameters characterizing the non-potentiality of the active region 
NOAA 11158 and its eruptive events using the magnetic field data from the He- 
lioseismic and Magnetic Imager (HMI) and multi-wavelength observations from 
the Atmospheric Imaging Assembly (AIA) on board Solar Dynamics Observa- 
tory (SDO). From the evolutionary study of HMI intensity and AIA channels, it 
is observed that the AR consists of two major rotating sunspots one connected 
to flare-prone region and another with CME. The constructed space-time inten- 
sity maps reveal that the sunspots exhibited peak rotation rates coinciding with 
the occurrence of the major eruptive events. Further, temporal profiles of twist 
parameters, viz., average shear angle, Oav, ttbcst, derived from HMI vector mag- 
netograms and the rate of helicity injection, obtained from the horizontal flux 
motions of HMI line-of-sight magnetograms, corresponded well with the rota- 
tional profile of the sunspot in CME-prone region, giving predominant evidence 
of rotational motion to cause magnetic non-potentiality. Moreover, mean value 
of free-energy from the Virial theorem calculated at the photospheric level shows 
clear step down decrease at the on set time of the flares revealing unambigu- 
ous evidence of energy release, intermittently that is stored by flux emergence 
and/or motions in pre-flare phases. Additionally, distribution of helicity injec- 
tion is homogeneous in CME prone region while it is not and often changes sign 
in flare-prone region. This study provides clear picture that both proper and 
rotational motions of the observed fluxes played significant role to enhance the 
magnetic non-potentiality of the AR by injecting helicity, twisting the magnetic 
fields thereby increasing the free energy, leading to favorable conditions for the 
observed transient activity. 
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Introduction 



Activity level of the Sun is defined by eruptive events like coronal mass ejections (CMEs) 
and flares occurring in the ARs consisting of groups of sunspots having positive and negative 
magnetic polarities. During the evolution of ARs, the sunspots exhibit motions, both proper 
and rotational, that lead to the storage of energy in the AR by increasing its magnetic non- 
potentiality. The eventual release of stored energy then occurs in the form of the eruptive 
events, such as, flares and CMEs. An investigation of the characteristics of sunspots can 
thus help towards the understanding of transient solar activities. 

Rota tional motions of and around sunspots have been observed by several authors in 
the past (lEvershedl llQlOl : iBhatnagarl 119671 : iMcIntoshI Il98ll ). More recent studies have also 
reported sunspot rotations about umbral center as well as about another sunspot with the 
availability of high s patial and temporal resolution data from space borne as well as ground 
based observatories flBrown et al.ll2003l : IZhang et al.ll2007l : lYan et al.ll2009l ). Rotational mo- 
tions can shear and twist the overlying magnetic field structures thereby storing excess 
energy. Vertical current, current helicity and transport helicity are some of the quantitative 
parameters derived from observations that are related to the twisting of fieldlines anchored 
in the moving footpoints. CMEs that disturb the space environment are thought to occur 
due to an over-accum ulation of hel i city a nd its ejection from the corona to the interplane- 
tary space. Recently, IZhang et al.l ( l2008l ) reported that helicity injection as inferred from 
the rotational motion of sunspots can be comparable to that as deduced for the entire AR 
by local correlation tracking (LCT) method. They also found good spatial and temporal 
correlation of sunspot rotational motion with two homologous flares. 

The consequences of sunspot rotations and related changes, both in photosphere and 
corona, have been investigated by several w'orkers. Relationships between the sunspot rota - 



tion and coronal co nsequences (IBrown et al 



sunspots, helicities ( lYan et al. 



with abnormal rotation rates 



2003 



flare productivity (lYan et al.ll2008l: IZhang et al.ll2008l ). statistics of the direction of rotating 



2008 



Tian fc Alexandei]l2006l : iTian et al.ll2008r ) 



and references there in"), and the association of flares 



Hiremath et al.ll2005l ) are some studies in this direction. 



In this study, we report on the rotating sunspots observed during February 11-16, 2011, 
in AR NOAA 11158. We present the results related to the changes in various parameters, 
viz., vertical current, twist parameters, helicity injection, and the free energy of the AR 



- 3 - 



during its evolution. Further, we examine whether the rotating sunspots were associated to 
the eruptive flare and CME events of the AR. Our aim is towards finding association of the 
observed uncommon rotation of the sunspots of this AR in terms of changes in various phys- 
ical parameters. We have organized the paper as follows: Observational data and analysis 
procedure is described in Section [2] which is followed by the results with possible discussion 
in Section [31 Finally, the summary of the study is presented in Section HI 



2. Observational Data and Analysis Procedure 

The observational data used in this study are obtained from the Helioseismic and 
Magnetic Irnager (HMI; ISchou et al.ll2012l ) and the Atmospheric Imaging Assembly (AIA, 
Lemen et al.l 120121 ) on board Solar Dynamics Observatory (SDO). 



AIA observes the sun in ten different wavelengths in soft X-ray, EUV and UV regimes 
to study the processes occurring at various atmospheric layers i.e., chromosphere, transition 
region and corona, at a pixel size of 0.6" and 12 second cadence. In addition, we also obtained 
high r esolution G-band 4 300A and Ca ll H 3970A images from Solar Optical Telescope 
(SOT; iTsuneta et al.ll2008l ) on board Hinode for understanding the associated chromospheric 
structures. Data on the CME and flare initiation times were obtained from the website^, B 

HMI observes the full solar disk in the Fe I 6173A spectral line with a spatial resolution 
of 0.5" /pixel and provides four main types of data at a cadence of 12 minutes: dopplergrams 
(maps of solar surface velocity), continuum filtergrams (broad- wavelength maps of the solar 
photosphere), and both line-of-sight (LOS) and vector magnetograms. Filtergrams are ob- 
tained at six wavelength positions centered at the Fe I line to compute Stokes parameters 
I, Q, U, V. These are then reduced with the HMI science data processing pipeline to re- 
trieve th e vector magnetic fie ld using the Very Fast Inversion of the Stokes Vector algorithm 
(VFISV; iBorrero et al.ll201ll ) based on the Milne-Eddington atmospheric mod el. The inher- 



ent 180° azimut hal ambiguity is resolved using the minimum energy met hod fjMetcaU 



1994 



Leka et al.ll2009l) . Data processing is carried out by the HMI vector team flHoeksema fc et al. 
2OI2I ) and the products are made available to the solar community. 



We used the high resolution continuum images of NOAA 11158 obtained from the HMI 
for the rotation measurement of sunspots. In order to estimate the twisting in various 
sites of the AR, we used vector magnetograms projected and remapped to heliographic 



^ http : //spaceweather . gmu . edu/ seeds/ 
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coordinates. The intensity images obtained from AIA were used to examine the coronal 
changes associated to the photospheric motions. All images both from AIA and HMI taken 
at different times were aligned b y differentially rotat ing to an image at central meridian 



by non-linear remapping method ( Howard et al. 199Cl[). Fo r our analysis, we have used the 



standard SolarSoft routines (SSW: lFreeland fc Handvlll998h implemented in Interactive Data 
Language (IDL). 

In order to further characterize the association of rotational motion on non-potentiality 
and in turn on the observed transient activity o f the AR, we examiri e the evolution of 
magnetic fluxes and the derived parameters (see iLeka fc BarnesI (120031 ) for definitions of 
various non-potential parameters) therefrom. We evaluate the vertical currents and shear 
or twist parameters to estimate the extent of non-potentiality. The complexity developed 
through the flux motions by helicity injection rate and consequent magnetic energy storage 
are also analyzed. Procedures for these calculations are briefly outlined below: 



1. 



3. 



We evaluate the total absolute flux in sub-regions Rl and R2 by summing over all 
the pixels F = X] l-^^l- "^^^ vertical current is computed from the Ampere's law 
V X B = /ioJ, where the permeability of free space fiQ = Att x 10~^ Henry/meter. 
Using the z-component of this relation, one can deduce the total absolute vertical 
current as = ^ |(V x 'B)z/fio\- 

Magnetic Shear: Non-potential nature of magnetic field is characterized by magnetic 
shear defined as difference between directions of observed transverse field and potential 
transverse field (lAmbastha et al.lll993t I Wang et al.lll994j ) and is given by 



S = B-9 = B 



cos 



1 B Bp 

BB„ 



T. 



where -B = |B| is magnitude of observed magn etic field an d Bp = |Bp| is magnitude of 
potential field computed from fourier method ( lGarylll989l ). The weighted shear angle 
in a region of interest with N pixels can then be obtained as 



WSA 



(2) 



and will estimate the overall angle of departure from potential field configuration. 

Alpha best (ahcst)'- The force free parameter abest is used as a proxy to monitor the 
extent of twist of the magnetic fieldlines in an AR. For this parameter, the extrapo- 
lated force free fields best fit the observed transverse fields in the sense of a minimum 
least-squared difference on the photosphere. The procedure of the fitting method for 
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calculating ojbest (jPevtsov et al.lll994t iHagino fc Sakuraill2004l ) is as follows. For a given 
a and with observed -B^, we compute linear- force free field whose transverse component 
is Bt cai- We then search for the a that gives minimum residual {R{a)) of transverse 
component in a range of values through the following equation 



Ra 



Y.[Bt,obs{x,y) - fBt^cai{x,y,a)f 



Yl -^t,obs 

where / 



t,obs 



(3) 
(4) 



The error in the fitting of residual and a that is effected in estimating ttbcst is given by 



^O^best 



(5) 



obs 



4. Average Alpha (aav)'- Another proxy for representing the twist is obtained by taking 
the avera ge of a from the force- free a ssumption of photospheri c fields, i.e., VxB = aB. 
Following iPevtsov et al.l ( 119941 ) and iHagino fc Sakurail ( l2004l ) , it is given by: 



E Jz{x,y)sign[B;,{x,y)] 
E\B.\ 



(6) 



The error in a^v is deduced from the least squared regression in the plot of B^ and J^, 
and is given by: 

r 2 ^ T.[Jz{x,y) - ae,^B^{x,y)]'^/\B,{x,y)\ 

{N-l)j:\B.{x,y)\ ^' 

where N is the number of pixels where \Bt\ > 150G (with 150G assumed as the lower 
cut-off for the transverse field.) 

5. Helicity Injection: One uses the heli city iniection rate to quantify the twist due to 
photospheric flux motions. Following iPariat et al.l ( 120051 ): 



dH 

~dt 



-1 
2^ 




[(x 



X u 



U')]n 



S J S' 



(8) 



where u(u') is the foot-point velocity at position vector x(x'), and B^ is the normal 
component of magnetic field from observations. This equation shows that the helicity 
injection rate can be understood as the summation of rotation rates of all the pairs of 
elementary fluxes weighted with their magnetic flux. Integration of above equation in 
the observed time interval provides the accumulated helicity, Aif = J [dH/dt)^t that 
is injected due to foot-point shear motions. 
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6. Magnetic Free Energy: One can get the estimate of non - potentiality by the magnetic 



virial theorem (j Chandr asekharl Il961t iMolodenskyI Il974 
total force-free magnetic energy as 



Lowl Il982l ) , which gives the 



Em. = / -^dV = — {xB-, + yB)BJxdy 
Jz>o 8vr Att J,^o 



(9) 



where B^, By are the horizontal and Bz the vertical components of magnetic field B 
at position (x, y) on the solar surface z=0. The applicability of this equation to the 
realistic situations like solar photospheric boundary must be restricted by force-freeness 
of the field and flux balance conditions. Despite its applicability, it should have some 
useful information about energy storage and release. Since the potential field state is 
a minimum energy state, by subtracting the potential field energy from the force-free 
energy one gets the upper limit of the free-energy available in the AR to account for 
the flares and CMEs. 



Howe ver, the photospheric field is not strictly force-free (IMetcalf et al.l Il995l : iGary 
200ll ) and flux balance in the region of interest are not fully satisfied. Therefore, 
one obtains incorrect estimates of the free-energy. To derive the error estimate for 
the force-free energy measurement, we have adopted the Pseudo- Monte Carlo method 
jMetcalf et al.ll2005r i. which gives the mean and variance of energy estimate by changing 
the origin of coordinate system in all possible ways. It includes both the statistical 
error as well as the error due to the departure from the force-free state of the fields. 



We have examined the errors given with the vector field data after inversion of stokes 
profiles. The maximum error in transverse field Bt in the entire data set varied in the 
range of 20-35G . The total field strength was observed to vary by 50G over a day due to 
various factors (IHoeksema &: et al.ll2012l ). Therefore, we carried out the calculations with 
\Bt\ > 150G and \Bz\ > 15G and neglected the pixels having magnetic fields below these 
thresholds, simultaneously. 

As require by Equation [HI the horizontal motions of photospheric fiuxes using LOS 
magnetic field dat a obtained at a cadence of 12 minutes by the method of Affine Velocity 
Estimator (DAVE; ISchuckll2006l ). It is a local optical fiow method that determines the mass 
velocities within the windowed region. Further, it adopts an affine velocity profile specifying 
velocity field in the windowed region about a point and constrains that profile to satisfy the 
induction equation. The maximum rms velocities of fiux motions in the entire AR were found 
to be distributed in the range 0.6-0.9 kms^^. These velocity maps of mass motions revealed 
strong shear fiows that included both proper and rotational motions. With the velocities of 
these fiux motions, we computed the helicity fiux density maps of the whole AR at pixels 
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above lOG to reduce computation time, then helicity injection rate, and accumulated helicity 
using Equation [81 

3. Results and Discussion 

The AR NOAA 11158 was a newly emerging region that first appeared as small pores 
on the solar disk on February 11, 2011 at the heliographic location E33S19. Subsequently 
it grew rapidly by merging of small pores to form bigger sunspots and developed to 
magnetic complexity on 2011 February 13. The HMl intensity maps in Figure [T] show the 
spatio-temporal evolution of AR 11158 during 2011 February 11-16. The prominent negative 
and positive polarity sunspots are labeled in frame (d) as SN1-SN3, SP1-SP3, respectively, 
along with the overlaid LOS magnetic field contours in frame (e) depicting the polarity 
distribution of the AR. The overall configuration of the AR is seen to be of quadrupolar 
nature. The proper motion of individual sunspots during February 13-16 are marked by 
arrowed curves traced by centroids of the sunspots as shown in frame (f). 

In what follows, we shall present evolution of the AR during 2011 February 13-15 using 
coronal observations in different wavelengths along with the changes in its magnetic and 
velocity field structures, the related physical parameters and the fiare/CME activities. 

3.1. Evolution and the Activity of the AR 

A careful examination of the animations of magnetogram and intensity maps revealed 
significant counter clock-wise (CCW) rotation of SNl and clock-wise (CW) rotation of SP2 
during February 13-15. We selected two sub-regions based on distinctive flux motions and 
associated transient activities. These are sub-region Rl covering SNl, and R2 around SN2, 
SN3, and SP2 as shown in the rectangular boxes in Figure (U^d). 

The spatial evolution in R2 displayed a large shearing motion of SP2 that itself rotated 
clock-wise during February 13-15. It then detached from SN2 and moved towards SP3 with 
small patches of both polarities appearing and disappearing over short periods of time. In 
Rl, a small positive polarity site SPl emerged to the north of SNl that rotated in a similar 
counter-clockwise direction along with a proper motion away from SNl. The rotation of SNl 
and its associated SPl twisted the field lines. These motions appears to have resulted in a 
large magnetic non-potentiality in this sub-region. This followed the first X-class event of 
the current solar cycle 24, an X2.2 fiare in R2 on February 15. In addition, several fiares of 
smaller magnitude, including four M-class and several C-class occurred during this period; all 
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associated to R2 except the M2.2 flare of 14 February which is observed to originate from Rl. 
Many intermittent CMEs were also launched from this AR. These flares and CMEs (listed in 
Table [T]) appear to be related to the observed motions in the AR. We have manually traced 
the the QuickLook AIA images to look for mass ejections in 304A channel locating the sub- 
region they belonged to and found that most of the CMEs originated from Rl (except the 
one associated with the X2.2 flare from R2). This inference is further confirmed by the 
STEREO space-craft informatioiifl. 



Table 1: List of Flares and CMEs 



AR Date 
(NOAA) dd/mm/yyyy 



Flares 

magnitude (time UT) 



CMEs 
(time UT) 



11158 13/02/2011 CI. 1(12:36), C4.7(13:44),M6.6(17:28) 
14/02/2011 C1.6(02:35),C8.3(04:29),C6.6(06:51) 
C1.8(08:38),C1. 7(11:51), C9.4(12:41) 
C7.0(13:47),M2.2(17:20),C6.6(19:23) 
C1.2(23:14),C2. 7(23:40) 
15/02/2011 C2.7(00:31),X2.2(01:44),C4.8(04:27) 
C1.0(10:02),C4.8(14:32), CI. 7(18:07) 
C6.6(19:30),C1.3(22:49) 



21:30,23:30 

02:40,07:00,12:50 

17:30,19:20 



00:40,02:00,03:00 
04:30,05:00,09:00 



Morphological changes observed in chromosphere, transition region, and coronal regions 
corresponding to some of the transient events are shown in Figure [2] using AIA observations 
in 304, 193 and 94A wavelengths. We further examine the characteristics of non-potentiality 
in the selected sub-regions. 

Figure [2] (al-a3) shows a large mass expulsion from sub-region Rl that turned into 
a large CME on February 14/18:00UT with the associated flare M2.2 at 17:20UT. It ap- 
peared to occur when the rotation of SNl had attained the maximum speed. This CME 
began at 17:30UT and continued well past 19:30UT as observed from the AIA/304A movies. 
The ejected plasma (marked by arrow in al) is visible in 304A as it is sensitive to the 
chromospheric temperature (~ lO^K). However, it is only partially visible in 193A and not 
discernible in 94A except for the cusp shaped fieldlines filled with the hot bright plasma at 
~ lO^K. Another such mass ejection of February 15/00:36UT that ensued from the same 
location is shown in Figure |2]^bl-b3). The LOS magnetic field contours are overlaid in frame 
b2 for reference, depicting the fieldline connectivities from the magnetic footpoints in the 



http : //spaceweather . gmu . edu/seeds/ 
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labeled sunspots. From the animations of magnetic and intensity images and coronal field- 
line connectivity, a significant role of new emerging flux at SPl can be inferred in these 
expulsions. 

Figure |2] (cl-c3) shows the intense X2.2 flare which occurred on February 15 in sub- 
region R2. This event was also associated with a halo-CME which began at 01:44UT. 
This flare produced abnormal mag netic polarity reversa l and Doppler velocity enhancement 
during the flare's impulsive phase (IMaurya et al.l 120121 ) . As mentioned earlier, the rotating 
flux footpoints, i.e., the sunspots, twisted the overlying magnetic structures in opposite 
directions. The rotational motion of the positive polarity spot SP2 about its negative polarity 
counterpart SN2 developed strongly twisted fieldlines along the polarity inversion line (PIL), 
clearly visible as a "sigmoidal" structure (c3). As reported ea rlier, such sigmoidal magnetic 
field structures are more likely to erupt (ICanfield et al.lll999l ). 



3.2. Measurement of Sunspot Rotation 

To de rive the rotat i onal p arameters of the sunspots, we followed the procedure as ex- 



plained in iBrown et al.l ( 120031 ). Essentially, we pursued the motion of a distinct feature in 
the penumbra as it is comparatively richer than the umbra. We uncurl the annular region 
of penumbra by transforming from the Cartesian (x-y) frame to the polar {t-6) plane. In 
this way, we can measure the position angle of the observed feature as displaced from its 
previous position as the sunspot rotated. We identified the centroid of the sunspot SNl by 
setting the 50% contrast between the umbra and penumbra in the time sequence of intensity 
images around the region of a circular disk of radius 18 arc-sec. This circular region was then 
remapped on to the polar coordinate system from the Cartesian, so that annular features 
appeared as straight strips from umbra to penumbra. Once this is done, any portion of 
the strip could be tracked in time to determine the angular distance by which a particular 
feature is displaced. 

We constructed the space-time diagrams for both the rotating sunspots, SNl and SP2, 
extracted at the radii of 11" and 7" from their umbral centroid, as shown in Fiugre [3](a-b), 
respectively. Rotational motion of the features in the penumbral regions can be seen as 
diagonal bright or dark streaks. The uncurling starts at westward point chord connecting 
the centroid position and the outer circle of annular portion and proceeds anti-clockwise 
about the sunspot centroid. These diagrams clearly show the counter-clock wise (CCW) 
rotation of SNl in sub-region Rl, and the clock wise (CW) rotation of SP2 in sub-region R2. 
A well observed streak is traced and plotted as marked by at the respective positions 
as angular displacement (6) and its time derivative i.e., angular velocity [dO/dt) with time 
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in Figure [3]^c and d). 

The sunspot SNl started rotating early on from February 14 which continued till Febru- 
ary 15/18:00UT. During the 24h period of February 14, it rotated by over 100° at an average 
rotation rate of 4°h~^ , peaking at 7°h~^ at around 23:50UT, i.e., just 3h before the X-flare. 
The timings of mass expulsions from this sub-region (cf.. Figure [2]^al)) also coincided well 
with the fast rotation rate of SNl. Except for the CME associated with the X2.2 flare, all 
other CMEs listed on February 15 were found to be associated with the sub-region Rl even 
in the decreasing phase of its rotation rate. This suggests that the new emerging flux SPl 
and large rotation rate of SNl might have played a role in the process of sigmoidal structure 
formation which gave rise to the CMEs that were often associated with flares. 

A similar analysis was carried out for evaluating the rotation of SP2 in order to examine 



its role in the X2.2 flare as reported by I Jiang et al.l ( 120121 ). For the sake of completeness and 



consistency, we have reproduced the results for SP2 in Figure [31 It exhibited rotation rates 
of 4.48°h~^ on February 14 before the X2.2 flare and 1.9 2°h~^ on February 15. These are 



consistent with the results deduced by iJiang et al.l (120121 ). Moreover, it rotated by 95° till 



February 15/04:00 UT which is also within the error limits of our manually tracked feature 
to their reported value of 107°. Notably, a small rotation rate of 0.9°h~^ was observed after 
the flare from the time profile showing that it rotated by a total of 25° till the end of the 
day. As is evident, there existed proper motion in addition to the rotational motion of SP2 
(Figure dl^f)). Therefore, it may not be appropriate to attribute the rotational motion alone 
to the X2.2 flare and other events observed with this region. 

It is thus evident that the two sub-regions consisting of the sunspots with large rotation 
rates were essentially the sites from where flares and CMEs originated in NOAA 11158. We 
speculate that the rotational motion led to the adequate storage of energy and injection of 
helicity which subsequently played the predominant role in the eruptions. In what follows, 
we shall discuss these aspects in further detail using the physical parameters deduced from 
the velocity and magnetic field measurements. 



3.3. Evolution of Physical parameters in sub-regions Rl and R2 

We now present a detailed description of the magnetic and velocity field observations and 
their derived parameters, as discussed in Section |21 in the sub-regions Rl and R2 containing 
the rotating sunspots. 

The vector magnetic field map of the sub-region Rl on February 14/18:00UT is shown 
in Figure IH^a) along with the map of tracked horizontal velocities of magnetic fluxes in 
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Figure Hl^b). The magnetic transverse vectors are seen nearly aligned with the polarity 
inversion line(PIL, thick dashed curve), implying a large shear along the SPl-SNl interface. 
The map of tracked velocities shows the velocity vectors exhibiting spiral or vortical patterns 
in the penumbral region of SNl as a fact of its rotation, and the velocity vectors aligned in 
the opposite direction to that of the magnetic field vectors. 

The new emerging positive flux SPl rotated in the CCW direction, along with a proper 
motion towards the main flare site R2. It is noteworthy that when the footpoints were 
dragged along one direction, the corresponding vector magnetic field pointed towards the 
opposite direction to retain its footpoint connectivity. In the process, the shear in magnetic 
structures built up. The shear motion along the polarity inversion line (PIL) is an effective 
mechanism for enhancing non-potentiality of magnetic field s. This can be identif ied as the 



alignment of transverse magnetic field vectors along the PIL (lAmbastha et al.lll993l ). Clearly, 
shear (Equation [1]) distribution map shows intense shear about the PIL with the rotation of 
SNl and SPl and proper motion of SNl as well. This process indeed enhanced the electric 
currents along PIL as is evident from th e the d istribution of calculated vertical current Jz in 



frame (c) of Figure HI Following IZhang] ( 1200 ll ). the Jz is decomposed to current of chirality 
Jf^ and current of heterogeneity J^, and examined their trend. The average ratio Jf^ / 
through out the time is deduced to be 1.11 implying the region is chirally dominated and 
is approximately consistent with force-free equilibrium. Interestingly, the shear current i.e. 

follows similar increasing trend with jf^ having correlation coefficient of 0.97, inferring 
that the rotational motion leads to increase both twist as well as shear contributing to total 
current Jz- Moreover, the force free parameter a = Jz/Bz is having negative distribution 
(of order 10^^ m~^) explaining the chirality with the counter rotation of the constituent 
sunspots in Rl. From the AIA observations of Rl, and also as confirmed by STEREO, it is 
observed that the intermittent mass expulsions turned into fast CMEs in this sub-region. 

A map of helicity flux density computed from horizontal flux motions is shown in Fig- 
ure |l](d). Through out the time interval, the helicity flux distribution is homogeneous with 
negative (dark) sign and implies left handed sense of chirality, further conforming to the 
observed CCW rotation of SPl and SNl. The bipolar flux system having such counter- 
clock (SPl) and counter-clock (SNl) rotations forms sigmoi dal patters which ar e efficient 



energy storage mechanisms to account for the flares/CMEs fICanfleld et al.lll999l ). This is 
also evidenced by the high resolution Hinode G-band continuum image showing fllamentary 
structures (frame (e)) that are consistent with the rotational motion of SPl and SNl. The 
corresponding Hinode chromospheric Ca II image (frame (f)) shows the flare ribbons along 
the PIL and overlying bright flare loops across the PIL (dashed line) associated with an 
M-class event of February 14. These observations suggest to infer that the rotational motion 
of SNl played a signiflcant role in creating favorable conditions for the eruptions. 
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The vector magnetic map of sub-region R2 is shown in Figure [5] at the time of X2.2 flare 
with a smoothed polarity inversion hne (thick dashed curve) separating positive and negative 
LOS magnetic field. This sub-region gave rise to several flares (listed in Table [T]) of varying 
magnitudes including the energetic X2.2 flare along with its associated CME. Importantly, a 
major portion of the flux (SP2 and SN2) is relatively associated to shearing motion, compared 
to flux within rotating SP2 itself, therefore this region R2 is shear motion dominated. As 
in Rl, the alignment of Bi vectors with PIL in this s ub-region als o showed strongly sheared 



magnetic fields. The tracked horizontal velocities (cf . . ISchuckll2006l ) displayed strong shearing 
motion of SP2 along PIL. These motions led to the stressing of field lines and possible flux 
cancellation as SP2 and SN2 possessed opposite fluxes. The distribution of vertical currents 
(Jz), another indicator of large non-potentiality of magnetic structures, too showed strong 
currents around the PIL in R2 in the form of J-shaped ribbons. High order shear is polarized 
along the PIL probably due to continuous shear motion of SP2 and rotation. This sub-region 
is roughly consistent with force-free equilibrium as the average ratio over the time Jf^j 
=1.29 indicating dominant chirality associated current. Nevertheless there exist dominant 
shear motion with magnetic shear and gradients, does not show increasing trend, but 
indeed there do observed in Jf^ as that in total current J^. It means that both rotational as 
well as proper motions are contributing to Jf^ the component that is parallel to to magnetic 
field. Therefore, the boundary motions in the evolution of the system that is in force-free 
equilibrium retains it to be in same force-free equilibrium. Distribution of helicity flux is 
inhomogeneous having mixed polarities and especially during the peak times of some flare 
events like M6.6 and X2.2, we noticed the negative helicity flux in the existing system of 
positive hehcity flux about the PIL. Injection of such opposite signed helicity is indicative 
of imminent transient activity. Whereas twisted penumbral flbril structure is seen in the 
Hinode G-band image on the photosphere, while, the bright chromospheric ribbons from 
X2.2 flare lying on either side of the PIL are seen in the Hinode Ca ll image. 

Temporal proflles of the deduced physical parameters for the sub-region Rl are plotted 
in Figure E] (left column). The total absolute flux doubled from 5 x lO^^Mx to 9 x lO^^Mx 
from February 13 to February 15. Thereafter, it showed no signiflcant increase. This increase 
of flux is gradual which is contributed by emergence of SPl and gradual evolution of SNl in 
a proportion of approximately 1:3. The absolute vertical current also increased along with 
the flux till February 15, and decreased afterwards. There was a pronounced increase in 
IJ^I, from 2.0 to 2.6xlO^^A, after 14/14:00UT. The average of shear (< S >) and weighted 
shear angle (WSA) are having similar trend, and explains the Jz temporal trend. This might 
be related to the increasing rotation rate of the sunspot, thereby increasing the shear of 
horizontal magnetic fields by dragging the foot-points accommodating large field gradients 
about the PIL. 
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The sign and magnitude of and abest are found to be consistent as both the param- 
eters are equivalent proxies representing the twist. The larger error bars in abest arise due 
to the 150G error in the transverse field given to propagate in the fitting procedure (See 
Eq. [7]). The negative sign indicates the left handed twist or chirality. This agrees well with 
the physically observed CCW rotation of the sub-region consisting of the sunspots SNl and 
SPl. The twist proxies, c^av and Obest, showed an increasing trend (in magnitude) with the 
rotation profile of SPl (see Figure [3]) further corroborating the increase in non-potentiality, 
as also evident from the absolute vertical current \Jz\- A similar trend also reflected from 
the temporal profile of helicity injection rate dH/dt calculated from the physically derived 
horizontal velocity field of the flux motions. The distribution of helicity flux is homogeneous 
with negative sign over this region Rl and the helicity injection rate increased (in magni- 
tude) with the rotation rate, reaching a maximum of —17.52 x 10^°Mx^h~^. Thereafter it 
decreased as did the rotation rate. As corona can not accommodate any amount of helicity 
that is being accumulated continuously, it would try to expel to still outer atmosphere in 
the form of CMEs so that the system to be stable, therefore explaining the CMEs observed 
from region Rl of this AR. A total helicity accumulation of —4.44 x lO'^^Mx^ was estimated 
in the sub-region Rl during this period. Since major part of the fiux(SNl) is associated to 
rotational motion, this accumulated helicity is mostly contributed by dominant rotational 
motion in this sub-region. 

The three parameters are different representations of twist or magnetic complexity ex- 
cept the manner in which they are estimated. Increased rotation rate of SNl seems to in- 
crease the twist so does the increasing trend of these parameters. The strong mass expulsions 
observed in AIA/304A at 12:50, 17:30, 19:20UT on February 14, turning into CMEs, fur- 
ther strengthen this evidence. This correspondence of the twist parameters Oav, «best,< S >, 
WSA, and dH/dt with the rotational profile of SNl suggests the predominant role of rotation 
in the observed expulsions/CMEs (refer to panels al and bl in Figure [2]). 

Thus we are able to identify the correspondence of the measured twist parameters with 
the observed rotation of the sunspot in sub-region Rl. However, the free-energy deduced 
for Rl did not show any obvious association with the observed CMEs. We believe that the 
broad peaks in free-energy profile seen during the CMEs (12:50, 17:30, 19:20UT) might have 
some relation but a rapid release of energy is not obvious. The maximum mean value of 
free-energy was estimated around 10"^^ ergs; adequate to account for the large eruptions. 

Temporal profiles of the deduced physical parameters for the sub-region R2 are plotted 
in Figure E] (right column) . Absolute flux for this sub- region was much larger in magnitude 
as compared to that in Rl. It increased monotonically from 10 x lO^^Mx to 14 x lO^^Mx. 
The magnitude of absolute current Jz too was correspondingly larger and showed appreciable 
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variations associated with local evolution. Peaks in the electric current profile are noticeable 
at the times of large flare events, i.e., M6.6 and X2.2. The currents decreased thereafter, 
implying the build up currents in the prefiare phase leads to relaxation of the magnetic 
field structure by releasing energy during flares. These peaks does not reflect in the < S > 
and WSA profile because of averaging property of these parameters. Note that WSA is 
higher by about 7° compared to that of Rl implying high level of non-potentiality over the 
region. Both the twist proxies, and Obest, agreed in sign and magnitude. Unlike in Rl, 
the positive sign of these parameters represents a positive or right handed twist in R2. A 
relatively large average value of at 8.5 x 10~^m~^ is suggestive of a stronger twist present 
in this sub-region. There seems no corresponding trend of any of these parameters with the 
rotational profile of SP2 probably complexity in the flux system is not only due to twisting 
but also shearing of fluxes. 

Unlike Rl, the sub-region R2 showed positive helicity injection at a much larger rate, 
with a maximum of 29.11 X lO^^Mx^h-^ The accumulated hehcity too was comparatively 
larger. A total helicity of 10.27 x lO^^Mx^ was injected during the 60h period by the observed 
flux motions and as a result the coronal helicity is likely to be positive. Within R2, shearing 
of flux system between SN2 and SP2 is dominated over twisting of single flux system of 
SP2, therefore it is likely that the accumulated helicity is mostly con t ribute d by the shear 
motion. This observance is consistent with reports of iLiu Sz Zhang] (j2006[ ) in AR 10488 
having shearing and twisting phases. They claimed that shear motion between two different 
flux systems indeed the major contributor of helicity injection in that AR. 

At the epochs of the M6.6 and X2.2 flares, sharp dips in helicity rate occurred cor- 
responding to the injection of opposite, negative helicity. Existence of such both signs of 
helicity flux in a s i ngle domain is supp osed to be sign of triggering transient flare events 
( iLinton et al.ll2001t iKusano et al.ll2004j ). However, such an observational evidence of nega- 
tive or opposite helicity flux may not be real change due to true helicity transfer, it could 
be flare related transient change effecting the magnetic fleld measurements. The observed 
distribution of negative helicity flux coincides spatially and temporally to the reports of pos - 
sible flare effects on Doppler and magnetic measurements in X2.2 flare flMaurya et al.ll2012l ). 
A speciflc study focusing the helicity inj ection and its behaviora l change during flares and 
CMEs from this AR 11158 in presented in lVemareddy et al.l (120 12[ ). They showed the change 
of helicity flux signal due to injection of opposite helicity could be a true change unless there 
is no impulsive effect of the of the flare on magnetic fleld measurements. In particular, the 
localized helicity flux distribution is affected by flare-related transient effects during these 
M6.6 and X2.2 flare cases having sudden dips in the helicity flux proflle as drawn in the 
figure. 
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The temporal profile of free-energy is plotted in the bottom frame of Figure |6l along 
with the GOES X-ray flux. The onset times of flares from the sub-region R2 are marked in 
this frame. The error bars represent standard deviation of free-energy obtained by changing 
the origin of coordinate system as explained before, representing the departure from force- 
freeness and flux balance. It is to note that the AR, as a whole, nearly fulfills the conditions 
for force-freeness and flux balance 0.05). Therefore, the constructed potential field from 
the extrapolation is not far from these conditions. The integration was carried out over the 
sub-regions of interest, i.e., Rl and R2, using eq O Although the force-free condition was 
nearly satisfied, a part of field lines from R2 were seen connected to Rl from the coronal AIA 
images. Hence, the condition for flux balance was not fully obeyed, posing some problem 
for the applicability of the virial-theorem. Therefore, the derived results are not conclusive 
enough unless the mean value of free-energy is much larger than the error bars. 

The most significant finding of the free-energy plot is the step-down decrease in free 
energy at the onset of the flares (indicated by the GOES X-ray profile). This provides an 
unambiguous evidence of energy release during the flares. Magnetic energy is believed to 
be stored by flux emergence and/or shear motions and is reflected in the form of electric 
currents. When an instability occurs in the field structure, the excess or free energy is 
released in the form of flares and CMEs. The flare events having such a correspondence are 
shown by vertical shaded bars. 

As the magnetic fluxes evolved and moved out of the sub-area R2, the fields departed 
increasingly from the flux-balance condition. Therefore, the free energy derived after the 
X2.2 flare are increasingly unmeaningful as evident from the large error bars. We estimated 
that the flux imbalance steadily increased from 5% on February 13 to 16% on February 
15. However, the mean value of free energy showed a clear step down decrease at the onset 
of the X2.2 flare, indicating an appreciable release of energy sufficient to account for this 
event. Therefore, the relatively large error bars notwithstanding, the variation of free-energy 
appears to be convincing at the time of this energetic X-class flare. 

The changes in potential and total energy estimated during the period from the onset to 
post-flare phase are given in Table [2] for some flares where we found the step-down decrease 
in the energy. The difference in available free-energy between the post-flare and the flare 
onset time giving the amount of energy released during the main flares, are described in the 
following. 

The X2.2 flare of February 15/01:44UT: The total energy Et ^ [4.92 ± 1.58] x lO^^ erg 
was estimated on February 15/01:48UT. It decreased to [3.67 ± 1.20] x 10^^ erg after the 
impulsive phase of the flare, i.e., at 02:12UT. The ratio Ex/Ep of the total energy with 
the energy of the respective potential state correspond to 2.01 and 1.46 at these two time 
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Table 2: Estimated total and potential energies during flares 



Flare Total Energy(i?r lO'^^ergs) Potential Energy(i?p lO^^ergs) 





Before Flare 


After Flare 


Before flare 


After flare 


M6.6 


1.11 ±0.39 


0.69 ±0.25 


0.61 ±0.22 


0.63 ±0.24 


X2.2 


4.92 ± 1.58 


3.67 ± 1.20 


2.43 ± 0.78 


2.47 ±0.79 


C8.3 


3.35 ±1.06 


3.16 ±1.00 


1.99 ±0.63 


1.97 ±0.63 


C6.6 


3.25 ±1.03 


3.20 ± 1.01 


2.02 ±0.64 


2.11 ±0.67 


C1.8 


3.32 ±1.05 


2.97 ±0.94 


2.11 ±0.67 


2.07 ±0.66 


C1.7 


2.27 ±0.79 


2.29 ±0.81 


2.88 ±0.82 


2.56 ±0.08 



instants. During this short 24 minute period, the excess energy released was estimated as 
[1.05 ± 1.04] X lO^^erg, suflicient to account for an X-class flare. 

The M6.6 flare of February 13/17:28UT: The ratio Et/Ep of 1.6 at the onset time of the 
flare at 17:24UT reduced to 1.08 at 18:12UT, i.e., after the flare. This implies the release of 
free-energy of [4.3 ± 1.7] x lO^^erg, i.e., consistent with the magnitude of the flare. 

Similarly, the free-energies released for the C-class events, listed in the Table [H shown 
by shaded bars and were estimated as (1.8 ±5.9) x lO^S (1.6 ±5.4) x 103\ (3.7 ±5.2) x 10^^ 
and (1.1 ±4.5) x lO^^erg. These estimated values are not statistically signiflcant as the mean 
values are smaller than the errors by a factor of two. Because the magnetic energy change 
is purely due to the change in the boundary conditions, it is expected that even minute 
changes in the magnetic fleld corresponding to a C-class flare would reflect in the mean 
values retaining error bars without much change. 

4. Summary and Conclusions 

This study suggests that the rotational motions of major sunspots in NOAA AR 11158 
were not only related to the transport of magnetic energy and complexity from the low 
atmosphere to the corona, but also appear to have played a major role in the onset of flares 
and CMEs. We have attempted to infer the relationship between the observed flux motions, 
variations in derived physical parameters in association with the dynamical transients, viz., 
flares and CMEs. A careful analysis of the flux motions leads us to better understand the 
dynamical nature of the photosphere at the interface between the convection zone where 
basic dynamo originates and the tenuous corona. 

The observations presented here provide a direct evidence for the energization of the so- 



-17- 



lar corona by the emergence and/or dynamics of magnetic flux-tubes tliat appear as rotating 
sunspots. This phenomenon is important for understanding how the solar atmosphere attains 
the conditions necessary for the release of energy and helicity observed in solar flares and 
CMEs. For a detailed study of various physical parameters characterizing non-potentiality, 
we selected two sub-regions, Rl and R2, of AR NOAA 11158, which consisted of rotating 
sunspots. In Rl with the rotating sunspot SNl, the estimated rotation rate is found to 
have a good correspondence with twist parameters ctav and ahest, < S >, and dH/dt. It 
showed that the intrinsic rotation of the sunspot increased the overall twist in sub-region 
Rl. Major CMEs in Rl occurred on 14 February at 12:50, 17:30, 19:20UT, which coincided 
well with the timings of large helicity injection rates and twist parameters. As the sunspot's 
rotation rate reached the maximum, the fluxes moved with maximum velocities to provide 
larger complexity in the connectivity. The calculated free-energies, however, did not show 
good correspondence with the observed CME events, except that broad peaks were observed 
during these CMEs. This is largely due to violation of fiux-balance condition(' 80%). It is to 



note that the occurrence of CMEs is also attributed to the new flux emergence (IMartin et al. 



19851 : IZhang et al.ll2008l ). However, the absolute flux profile of Rl does not show such sudden 



emergence of flux related to emergence of SPl. 

No clear relationship of the rotation rate was found with the twist parameters and 
dH/dt for the sunspot SP2 in R2. Shear motion of SP2 (cf.. Figure [T]f) seems to be the 
dominant factor, which obliterates the rotational effect of SP2 on the twist parameters and 
helicity rate. This shear motion caused the storage of magnetic energy by stressing the 
field lines. The intermittent release of energy during the observed flares is discernible. A 
total free-energy of (1.05 ± 1.04) x lO^^erg is released during the X2.2 flare. Th is is three 



times larger than the derived energy from coronal magnetic field extrapolations by lSun et al. 



(120121 ) over the entire AR. Spectral line-profile reversals occurred during the impulsive phase 



of this flare af f ecting the magnetic and velocity field measurements, as reported recently by 



Maurya et al.1 (120121 ). As a result, the injected helicity rate showed an impulsive negative 
peak by the sudden appearance of negative helicity density indicating flare-related transient 
effect (cf.. Figure [Sji). 

The M6.6 flare of February 13/17:28UT was another large but less impulsive event 
which released (4.3 ± 1.7) x lO^^erg of energy in the 48 minute period. Step down decrease 
of free-energy was found not only during the two large flares, but also in some flares of 
smaller magnitude. This reveals the storage and release of energy in the flares occurring in 
sub-region R2. It is to mention that we have not carried out the analysis of errors due to the 
spectro-polarimetric as well as random noise and their effects on the free-energy estimation, 
as there is no established method available on the data set as yet. Moreover, the estimation 
of free-energy suffers by the departure from the validity conditions of virial theorem due to 
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the selection of small sized areas of the sub-regions. 

The accumulated helicity contributed by the rotational motion of the sunspot can also 
be obtained as 

AH = Ae^\ (10) 

271 

where $ is the magnetic flux and A6 is the total rotation angle of the sunspot. Since 
the sunspot SNl rotated by ~ 160° during the 60h period of our study, and assuming the 
average magnetic flux of this sunspot as 3.5 x lO^^Mx, the accumulated helicity by rotation 
is estimated as -5.44 x lO^^Mx^. This is consistent with the derived value of -4.44 x lO^^Mx^ 
from the tracked velocities of the fluxes. 

Thus, this study demonstrates that both shear (dominated in R2) and rotational (dom- 
inated in Rl) motions of the observed fluxes enhanced the magnetic non-potentiality of the 
active region by injecting helicity. This helped in increasing the free energy of the AR's mag- 
netic field by increasing the overall complexity leading to the conditions favorable for the 
eruptions. The deduced physical parameters describing the level of magnetic non-potentiality 
and the free energy showed a reasonably good correspondence with the observed transient 
activity of the AR. This study also provides a clear signatures of energy release during the 
major, energetic transients of the AR. 

We thank an anonymous referee for the comments that helped us to improve the read- 
ability of the manuscript. The data have been used here courtesy of NASA/SDO and HMI 
science team. We thank the HMI team for making available the processed vector magnetic 
field data. The authors are grateful to Dr. Xudong Sun of Stanford University for his help 
in handling the vector magnetogram data. This work utilizes the Ca ll H and Continuum 
data from the Solar Optical Telescope (SOT) on board Hinode. 
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Fig. 1. — HMI intensity images showing the evolution of AR NOAA 11158 during six day 
time period. The main sunspots are labeled by SP/N* in (d) along with the LOS magnetic 
field contours overlaid in red (blue) color at 150(-150)G levels in panel(also in subsequent 
figures unless specified) (e). Proper motions of individual sunspots are traced along the 
arrowed curves as in panel (f). The two rectangular boxes mark the selected sub-regions Rl 
and R2 for further study. 
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Fig. 2. — The coronal activity observed during the evolution of the AR in AIA wavelengths. 
(al-a3): A mass expulsion (dashed arrow) from Rl which eventually turned to a CME on 
2011 Feb 14/18:00UT. (bl-b3): Another mass expulsion observed on 2011 Feb 15/00:36 UT 
from the same location, (cl): Bright emission seen during the X2.2 flare, saturating the 
I93A detector (c2) and the twisted flux ropes along the polarity inversion line as seen in 
94A (c3). 
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Fig. 3. — Space-time map of a strip in (a) the sunspot SNl (at a distance of 7 arc-sec 
from the sunspot's centroid), (b) sunspot SP2 (at a distance of 11 arc-sec from the sunspot's 
centroid), after remapping the annular region of the sunspot to the t-6 plane. Rotation angle 
{6) and rotation rate (dO/dt) of the sunspot's umbrae are inferred from the well-appearing 
diagonal feature in (a)-(b) sampled along the path marked by and plotted with time 

in (c)-(d), respectively. 
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(g) Hinode G band 



(h) Hinode Ca II 



Fig. 4. — Sub-region Rl on February 14/18:00UT: (a) transverse magnetic field (Bf) vec- 
tors overlaid on LOS magnetic field drawn with polarity inversion line (dashed curve), (b) 
horizontal velocity vectors showing vortical-like pattern due to rotation of SNl, (c) verti- 
cal current (J^ = (V x B)z/fio) distribution, (d) shear angle map indicates strong shear 
about the PIL of SPl and SNl with their rotation and shear motion of SPl, (e) Alpha 
(a = (V X 'B)z/Bz) distribution with iso-contours of LOS field indicates overal negative 
twist of the region (f) helicity flux density map showing negative helicity distribution in the 
sunspots that is consistent with negative "alpha" , (g) Hinode continuum G-band image, and 
(h) Hinode chromospheric Ca ll image with the bright flare ribbons of M2.2 at 14/17 : 20UT 
around the PIL. 
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Fig. 5. — Sub-region R2 on February 15/02:00UT: (a) transverse magnetic field {Bt) vectors 
overlaid on LOS magnetic field with PIL separating positive and negative polarities, (b) 
tracked horizontal velocity vectors overplotted on LOS magnetic field, (c) Vertical current 
distribution ( J^) showing intense currents accumulated on either side of PIL with the shear 
motion of SP2, (d) Shear angle map also implies strong shear polarized about the PIL with 
the continuous shear motion of SP2, (e) Distribution of twist parameter "alpha" indicates 
positive twist of the region, (f) helicity flux density map with intense negative helicity flux 
about PIL probably due to fiare-effects on magnetic field measurements, (g) Hinode contin- 
uum G-band, and (h) Hinode chromospheric Ca ll image showing the flare ribbons of X2.2 
around the PIL (dashed curve). 
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Fig. 6. — Temporal evolution of various physical parameters characterizing the non- 
potentiality in sub-regions Rl (left column) and R2 (right column). Arrows in the lower- most 
panel of the left column mark the CME timings while those in the right column correspond 
to the GOES-15 flares. The three major flares are marked by dashed vertical lines in all 
panels on the right column, for reference, while the shaded vertical bars indicate the release 
of free-energy in step with the onset of flares. Note that all the twist parameters Oav, «bcst, 
average shear (< S >), and dH/dt has clear correspondence with the rotational profile of the 
sunspot SNl implying that sunspot rotation has direct role in increasing the non-potentiality 
in sub-region Rl. 



